Nonlocality and entanglement in the XY model 
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Nonlocality and quantum entanglement constitute two special features of quantum systems of 
paramount importance in quantum information theory (QIT). Essentially regarded as identical or 
equivalent for many years, they constitute different concepts. Describing nonlocality by means of 
the maximal violation of two Bell inequalities, we study both entanglement and nonlocality for two 
and three spins in the XY model. Our results shed a new light into the description of nonlocality and 
the possible information-theoretic task limitations of entanglement in an infinite quantum system. 
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Schrodinger 's reply pQ to the paradox posed by Ein- 
stein, Podolsky and Rosen (EPR) [? motivated the mod- 
ern notion of entanglement in a quantum system. EPR 
suggested a description of nature, called "local realism" , 
which assigned independent properties to distant parties 
of a composite physical system, to conclude that QM 
was an incomplete theory. Schrodinger, instead, did not 
recognise such conflict and regarded entanglement as the 
characteristic feature of QM. 

The most significant progress toward the resolution of 
the EPR debate was made by Bell [3]. Bell showed that 
local realism, in the form of local variable models (LVM), 
implied constraints on the predictions of spin correla- 
tions, known as Bell inequalities. Spatially separated ob- 
servers sharing an entangled state and performing mea- 
surements on them may induce (nonlocal) correlations 
which cannot be simulated by local means (violate Bell 
inequalities) . This limitation to our physical understand- 
ing is nowadays exploited for implementing information- 
theoretic tasks. 

Ever since Bell's contribution, entanglement and non- 
locality were essentially identified as the same thing. 
With the advent QIT, interest in entanglement dramat- 
ically increased over the years for it lies at the basis of 
several important processes and applications which pos- 
sess no classical counterpart [H [S] . 

Confusion between nonlocality and entanglement arose 
when the usefulness of quantum correlations was put in 
doubt (see [O]). The nonlocal character of entangled 
states was clear for pure states since all entangled pure 
states of two qubits violate the CHSH inequality and are 
therefore nonlocal {Gisin's theorem) [7]. However, the 
situation became more involved when Werner [S] discov- 
ered that while entanglement is necessary for a state to 
be nonlocal, for mixed states is not sufficient. 

Entanglement is commonly viewed as a useful resource 
for various information-processing tasks. Yet, there exist 
certain procedures, such as device- independent quantum 
key distribution [2j and quantum communication com- 
plexity problems |10j , which can only be carried out pro- 
vided the corresponding entangled states exhibit nonlo- 
cal correlations. Therefore we are naturally led to the 



question whether nonlocality and entanglement consti- 
tute two different resources. 

The purpose of the present work is to shed some light 
upon the relation between entanglement and nonlocality, 
through the maximal violation of a Bell inequality, in an 
infinite system, namely, the XY model [TT] . 

The model Hamiltonian of the anisotropic one- 
dimensional spin-i XY model in a transverse magnetic 
field h for N particles is given by 
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are the Pauli spin-^ 



where ct^ — 25^ {u — x,y,z) are ttie Fauii spm-| op- 
erators on site j, 7 G [0,1] and (T^"*"^ — al^. The 
XY model (IT]) for iV = cx) is completely solved by ap- 
plying a Jordan- Wigner transformation [TT| [T2j, which 
maps the Pauli (spin 1/2) algebra into canonical (spin- 
less) fermions. This model undergoes a paramagnetic- 
to-ferromagnetic quantum phase transition (QPT) [13] 
driven by the parameter h at he = 1 and T = 0. 

We will provide evidence for an anomaly that regards 
entanglement and nonlocality in the XY model. To such 
end, we will consider the correlations existing between 
two sites or qubits (bipartite case) and three sites or 
qubits (tripartite case). 

Two qubits The general two-site density matrix is ex- 
pressed as 



{R) 1 



/ J uv u 



(2) 



R — j — i indicates the distance between spins, 

{ R) 

{u,v} denote any index of {ctqi ctx, Cy, CTz}, and Tuv = 
(crjj (81 al). Due to symmetry considerations, only 
{rif ^ , T^y'' , rif ^ , ri^^ } do not vanish. Barouch et al [H] 
provided exact expressions for two-point correlations, to- 
gether with all the dynamics associated with an external 
h{t). Let us consider the case where h jumps from and 
initial value ho to a final value ft,/ at t = (the equi- 
librium case is easily recovered when hf = /lo) and the 



R — 1 configuration. Following [T5], one obtains that 
Gri, Izz ^= ^0 — CziGr—i — oiS-^i and 
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Si, where 

7 r tanh [i/3A(feo)] 

: - / d(piiu\{R(p) N.,^, — ^ 

7r7o A(;io)A2(/i/) 

X [7^ sin^ (/) + (/iQ — cos 0) [hf — cos 0) 

~(ft.o - hf){hf - cos 0) cos(2A(/i/)i)] 

1 r ,, ,_,tanh||^A(Ml 

71" Jo A('io)A^(/i/) 

X [{7^sin^ 0+ (ft-o — cos 0) (ft./ — cos 0)} (cos 0— 
ft/) - (fto - ft/)7^sin2 0cos(2A(ft/)t)], (3) 



5p = 



7(fto - ft/) 



■ /n^N • j.sin[2A(ft/)t] ,^, 
-"'^(^^)"^^^ A(MA(ft/) ' (^) 




with A(ft) = [7^ sin^ 4> + {h - cos f/")^]^/^. Gr is the two- 
point correlator appearing in the Wick theorem calcula- 
tions, and Mz = ^Gq- 

Most of our knowledge on Bell inequalities and their 
quantum mechanical violation is based on the CHSH in- 
equality [7] . With two dichotomic observables per party, 
it is the simplest nontrivial Bell inequality for the bi- 
partite case with binary inputs and outcomes. Quan- 
tum mechanically, these observables reduce to Aj(Bj) = 
aj(bj) • a, where aj(bj) are unit vectors in M'^ and a — 
(cxjtTyjCTz) the Pauli matrices. Violation of CHSH in- 
equality requires Tr{pi- Bchsh), that is, the expecta- 
tion value of the operator Bchsh 
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FIG. 1. (colour online). The oscillating curve depicts Bq'^sm 
([7| versus time (in seconds) for /iq ~ 0.5, hf — and 7 = 0.5 
for the evolved state p\- (t) |2| at T=0. Final and equilib- 
rium values do not coincide (non-ergodic). Inset a) depicts 
BcfFsM versus h for 7 = 0.5 and different distances between 
spins R — f , 2, 3 and 00 (dotted line). Though no Bell viola- 
tion is observed, B^flgj^f exhibits a long range order. Inset 
b) depicts equilibrium nonlocality (upper curves) and con- 
currence (lower curves) measures versus h for 7 G [0, 1] . The 
isotropic case 7 = collapses to 2 for ft, > 1 . See text for 
details. 
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Ai(8)Bi-|-Ai®B2+A2®Bi-A2«)B2 (5) where T^ = 



(rpiR) rr.{R) ^(fl.) ^ , ■ 



to be greater than 2. We shall take the optimum value 
over all {aj,bj} as a proper measure for nonlocality in 
our state p^ ' (|2|. This procedure is presented in de- 
tail elsewhere [Hj- Given a general two qubit state p in 
the usual computational basis, we change it into the well 
known BeU basis {1$+), |$"), |*+), |*")}. The ensuing 
matrix p = pu +p± is decomposed into two contributions, 
where only terms in pu 



mm{[T!,i^'']\[T^'^']\m^>Y). Fig. 1 depicts the evo- 
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contribute to Tr(pBcHSH)- In the XY model, state ^ 
is almost Bell-diagonal except for pfj — 5^^2/1 which is 
null in equilibrium (ft/ = ho). 

Optimization is carried out and after some algebra, we 
obtain 2V^y/{pii - PaaY + (P22 - ^33)^ + '^{p[2Y, with 
the diagonal elements of ^ arranged so that pn > P22 > 
P33 > P44- The specific form for our state ^ reads as 

B^SlfsM ^ max Tr{pf Bchsh) 



lution of nonlocality (|7|) for the case 7 — 0.5 and (fto = 
0.5, ft/ = 0). For 7 = 1 (Ising) and (fto = ftc,ft/ = 0) 
([7]) would oscillate indefinitely. This nonlocality measure 
oscillates around a value distinct from the one expected. 
Just as in the case of the Mz in [12] , or entanglement in 
the XY model [15] (as measured by the concurrence C), 
([7]) does not reach its equilibrium value, which entails 
that nonlocality is also a non-ergodic quantity. 

The equilibrium case (ft/ = fto) is considered in Fig. 
la. Bc'hsm evolves for different configurations up to 
R — 00 [12]. It is apparent that nonlocality exhibits a 
long range behavior, while C rapidly tends to zero [16j . 
Comparison with entanglement appears in Fig. lb for 
all 7-anisotropies. The usefulness of (bi)entanglement 
between spins in the XY model is questioned in quantum 
information processing by the fact that the concomitant 
correlations never violate a Bell inequality (B™|fgJ^^ < 
2Vft,7). 

All previous quantities are ultimately described in 
terms of several Gr, so that they all diverge at ft = 1 
in the same manner. Let us consider for simplicity 
Mz{h) = iGo = |^5LJ-d</.[72sin2<^ + (ft-cos</.)2]i/2. 

For 7 = 1 we have M,(ft) = |^(?(^i?[|^]) = 



^ ['^K{j^^) + '^E{j^^)] , where KiE) is the com- 
plete eUiptic mtegral of the first (second) kind. jf^Mz 
diverges logarithmically at h = 1 following the diver- 
gence of K, and so does C and B^lfgj^j, including non- 
equilibrium [t = oo) values. Therefore entanglement and 
nonlocality exhibit one additional feature beside non- 
ergodicity: they both signal the presence of a QPT at 
T=0. 

Three qubits Nonlocality in the three qubit case is ex- 
plored through the violation of the Mermin inequality 
[17j . The Mermin inequality reads as rr(/9i?Mermin) < 2, 
where BMermin is the Mermin operator 
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Mermin 



B„ 



Baib2b3 - -Bbiaaba ^ -Bb^fe^Qj, (8) 



with B„ 



= u-(T®v-(T(8)w-(7 with a = {ox-, Oy^Oz^ 
being the usual Pauli matrices, and aj and bj unit vectors 
in M'^. Notice that GHZ states maximally violate the 
Mermin inequality. As usual, we will employ 



Mermin" 



max Tr{pB Mermin) 
ai,b; 



(9) 



as a measure for the nonlocality of the state p. 

The XY model is completely solvable, a fact that allows 
us to compute - as in the previous case of two sites - 
the reduced density matrix for three spins without the 
explicit construction of the global infinite state of the 
system. The reduced state of three spins reads as 






Erp(a,\ 
uvu 



(10) 



where i < j < k indicate the positions of the three spins 
and a = j — i,b = k — j their relative distances, {u, v, w} 
denote indexes of the Pauli matrices {cro,ax,cry,az}-i and 
^i®'^w)ab- Similarly to the calculation of 



rp(a,b) _ I I 



the two-spin correlations computed by Barouch et al [H] , 
based in turn on the work by Lieb et al |11| . we extend 
them to the three party case by using the well known 
Wick theorem in quantum field theory. Due to the sym- 
metry of the XY model, some of them vanish. Further- 
more, as far as nonlocality is concerned, among those cor- 
relations who survive only four of them contribute to d9| , 
namely, Ttz\Ttx\Ttx^ and Tt^\ These three-spin 
correlation functions Ti^x , Tx%'z and Tiiz are given by 
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FIG. 2. (colour online). Fig. a) depicts the value of 
A/ermm™^" versus h ([9| for several 7 values for the a = 
6=1 block configuration. Fig. b) shows similar curves for 
a = l,b = 2 and a = 2,b — 2, with 7 — 0.5. Fig. c) depicts 
the evolution of multipartite entanglement Svn {a = b = 1) 
versus h for several values of the anisotropy. Fig. d) plots 
Mermin'^'^^ versus SvN for 7 = 0.5 and all previous config- 
urations a = fe=l, a = 1,6 = 2 and a — b = 2. Notice the 
apparent monotonic decreasing evolution. See text for details. 
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symmetry. 

Correlations between parties strongly depend on their 
relative positions a and 6. We shall distinguish two types: 
i) the one forming a block of three in consecutive sites 
(a = 6 = 1), as well as a centered configuration [a = h — 
2), and ii) two nearest-neighbors spins plus an additional 
one (a = 1, 5 = 2). 

Optimization of M ermin™^'^ ([9]) for any configuration 
of the spins along the chain is carried out as follows. 
Once the observers' settings {aj, bj} are parameterized in 
spherical coordinates (sin^/j cos0fc,sin0fc sin0fe,cos0fc), 
the problem consists in finding the supremum of Q over 
the set of {fc = 1..12} possible angles. Here, too, we will 
consider ground state nonlocality. 

Let us consider the block configuration a = h = \. 
The orientation of the settings {aj,bj} that max- 
imizes Mermin'"'^^ ^ is such that {x — z plane 
(11) only, that is, we deal we real qubits for u^ = OVu) 
{as = ai,bi = a2,b2 = -ai,b3 = a2} for 
the whole range of h. Explicitly, Mermin'"^^ = 
[ai((ai)2 - 3(af)2)rili^) + (ai((af)2 - 
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ensuing analytic form of Mermin^^^ versus h - a 
complicated rational function with radicals- is cal- 
culated by recourse to convex optimization. For the 
sake of generality, let us also consider the centered 
configuration a = b = 2. After some algebra, we obtain 
that one disposition of the observers that provides an 
analytical expression for a lower bound to ([9]) is given by 
(again in the x — z plane) {ai — (sin0aj,cos0Qj^),a2 = 
(0,-1), ag = (sin6'ai,- cos 6*01 ),bi = a3,b2 = 
(1,0), bs = — ai}- Hence, we obtain that 

max0^^[2cos'^ea^T^fz^ + 2 sm0a^ cos 9 a^TJ^li^ 
2 sin" gg.rJL'^ + 2sin^,^ cos Oq^TJl'^] < Mermin '"'^^ < 

^^{Ttl^f + ^T^hy + A{Tt.y + 4(riL"V. The 
analytic form of the lower bound - not provided here - 
is of the same nature of that of the a = b = 1 case. In 
point of fact, the lower bound becomes an equality for 
all {a,b) shortly before h — I. Additionally, the upper 
bound also applies to all configurations. 

Fig. 2a, 2b present the situation, where Mermin^^^ 
is depicted for different values of the 7-anisotropy and 
different configurations of the spins. Numerical calcu- 
lations perfectly agree with the corresponding analytic 
expressions. As h grows, the state pfA approaches 
I 4-4-1- ) ( 4-4-4- I as expected (ferromagnetic phase) , and 
never violates the Mermin inequality, which entails an 
inherent limitation to the usefulness of entanglement it- 
self. 

Characterization of entanglement is of paramount rel- 
evance in QIT |18| , yet no operational measure is avail- 
able to date that quantifies genuine multipartite entan- 
glement. We will nevertheless employ the sum of the 
von Neumann entropy of the reduced density matri- 
ces of the three spins of py^f, , SyN — 3S'.ujv(pi), with 
p. = 1(1+ {(Tz)fTl), {(Tz) =Go = 2Mz. We shall consider 
ground state entanglement (T=0), as well. The specific 
form of SyN is depicted in Fig. 2c for several values of 
the 7-anisotropy as a function of h. The monotonic de- 
creasing tendency of entanglement is apparent for any 
7- value. As h grows, the fidelity between p^-J^. steadily 
tends to 1. 

Remarkably, we observe opposite tendencies revealed 
by entanglement S^m and nonlocality ^. Recall that our 
model undergoes a second order QPT in the ground en- 
ergy. In point of fact, since jf^SyN oc jjiMz, this measure 
of multipartite entanglement (its derivative) diverges log- 
arithmically at he — 1. Surprisingly, our measure of 
nonlocality (|9| also displays such divergence, along with 
its bounds. In Fig. 2d we show the dependency of 
Mermin"^'^^ versus entanglement SyN for several values 
of 7 in the (a = 1,& = 1) configuration. Finally, we en- 
counter for the first time a multipartite system where not 
only is entanglement (its first derivative) a good indicator 
of a QPT, but also nonlocality. 



Conclusions We have studied how nonlocality - 
measured by the maximal violation of a Bell inequality- 
compares to entanglement in a condensed matter system. 
Although two instances (two and three sites) have been 
considered, our results may properly generalize to any 
block of spins. For the bipartite case, we have computed 
the exact value of B^ff^j^j during time evolution and 
in equilibrium. In either cases our nonlocality measure 
([7| displays a non-ergodic behavior, is able to detect a 
QPT, and limits the QIT-related tasks involving bipar- 
tite entanglement along the infinite chain since no Bell 
inequality is violated. A similar situation occurs in the 
tripartite case, where non-violation of local realism in 
the XY model takes place as well. Also, entanglement 
and nonlocality both indicate a QPT yet exhibit oppo- 
site evolutions in the phase diagram. Finally, nonlocality 
can also constitute a complementary resource in infinite 
quantum systems. 
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